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Abstract 

We modify the standard Wilson action of 517(3) lattice gauge theory by adding an extra 
term which suppresses colour magnetic currents. We present numerical results of simulations 
at zero and finite temperature and show that colour magnetic currents strongly influence the 
confining properties of 517(3) lattice gauge theory 


1 Introduction 

The mechanism which leads to permanent confinement of quarks at low temperatures is an 
interesting and still open question in gauge theories. A promising conjecture which could explain 
the existence of a non-zero string tension is the hypothesis that the QCD vacuum behaves 
dually to a superconductor [4j] . In this model colour magnetically charged particles condense in 
the ground state of QCD and squeeze the colour electric field between a pair of static quarks 
into a narrow flux tube leading to a confining potential. A possible realization of the dual 
superconductor picture in non-Abelian SU(N ) gauge theories was proposed by’t Hooft in the 
framework of Abelian projection |2j]. In this approach, the physically relevant degrees of freedom 
are identified by fixing the gauge up to the largest Abelian subgroup U( l)^ -1 of the SU(N) 
gauge group. After gauge fixing one is left with an effective Abelian theory which contains in 
addition to photons magnetic monopoles. The monopoles arise as singularities in the gauge fixing 
condition and their condensation is expected to be crucial for the existence of a non-zero string 
tension. There is a lot of support for the Abelian projection approach from numerical results 
of lattice simulations Q. For instance, the Abelian gauge fields obtained after gauge fixing 
carry almost the whole asymptotic string tension Q. This observation is referred to as Abelian 
dominance. However, most of the numerical results are dependent on the choice of gauge which 
fixes the non-Abelian degrees of freedom ||], and a dual superconductor is convincingly observed 
only in one particular gauge, the so-called maximal Abelian gauge j|]. The question arises 
whether this apparent ambiguity in the gauge fixing procedure can be avoided by investigating 
the dual superconductor picture in terms of the original non-Abelian degrees of freedom. 

A first attempt in this direction was made in ref. [!(| where dual superconductivity in SU (3) 
lattice gauge theory was studied numerically in a gauge invariant formulation. An operator for 
the magnetic current was introduced allowing the determination of its curl in the presence of a 
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static quark antiquark pair without need of fixing the gauge. At finite temperature the result of 
a numerical simulation is the following Q : In the confined phase the curl of the magnetic current 
is proportional to the electric field indicating the validity of a dual London relation. This seems 
to be in agreement with the dual superconductor picture of confinement. In the deconfined phase 
the curl of the magnetic current vanishes and the electric field shows Coulomb-like behaviour. 
At least qualitatively the results presented in ||6!] agree with the results reported in ref. @] which 
were obtained in the maximal Abelian gauge. 

In this letter, we continue the gauge invariant investigation of the dual superconductor 
picture of confinement and study the role of the magnetic current operator introduced in ref. || 
in the pure gluonic vacuum. We introduce a chemical potential A and add an extra term to 
the standard gauge field action suppressing large magnetic currents. By calculating the string 
tension and the critical coupling for the finite temperature phase transition we will show that a 
suppression of magnetic currents strongly influences the confining properties of the considered 
SU( 3) theory. 

2 The Model 

We consider pure SU (3) lattice gauge theory on a four-dimensional Euclidean lattice of spacing 
a with periodic boundary conditions in space and time direction. An action convenient for many 
purposes is the standard Wilson action Q which for the SU (3) gauge group reads 

Sw = P E (l —^ReTrCV(x)) , (1) 

where /? = 6/g 2 is the inverse coupling and U pu (x) is the product of link variables U p £ SU(3 ) 
around an elementary plaquette in /ii/-direction at lattice site x: 

U^ v {x) = U^(x)U v (x + (j)Ul(x + v)Ul(x). (2) 

According to a derivation of Gauss’ law on the lattice |j] the field strength F pv (x) corresponding 
to the action ([l]) is given by the expression 

ga 2 F*„(x) = F*„{x)= ^ (u pu (x) - U^xfj , (tl) = traceless. (3) 

The upper index x in the notation of F*(x) states the local colour coordinate system in which 
the field strength is measured. It is changed to a system defined at an arbitrary lattice site x' 
by performing a parallel transport of the field strength along a Schwinger line connecting the 
site x' with the original site x. 

Let us now introduce the quantity Q given by the expression 

(4) 

where we use the discretized covariant derivative of the field strength 

A+ 0) - J%(x) = U u (x)F^(x + v)Ut(x) - F* pa (x). (5) 

On the lattice definition ([!) has a simple geometrical interpretation. Each component of J mtfl 
corresponds to a three-dimensional cube built of six plaquettes measuring the flux out of the 
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cube. We therefore call J mtfl colour magnetic current. The covariant derivative in Q) guarantees 
all contributions to the total flux out of the cube to be taken in the colour coordinate system 
at x. As can be easily shown || the current (|j) fulfills a conservation law in terms of covariant 
derivatives. 

In ref. Q we already discussed the definition of the quantity J m /t concerning the validity 
of the Bianchi identity. We think it is necessary to extend this discussion and add at this 
point some further comments concerning especially the general form of the Bianchi identity in 
a lattice formulation of SU( 3) gauge theory. In a continuum formulation the right hand side 
of ([|) corresponds to the Bianchi identity and the current J m ,n vanishes. On a lattice with 
finite lattice spacing a, however, definition (^) may not be identified with a lattice version of 
the Bianchi identity. This can be understood as follows: As was shown in ref. [;|] the partition 
function Z of SU( 3) lattice gauge theory expressed in terms of link variables U fl (x) £ SU{ 3) 

Z = J V[U^x)\eMSw[U^}) ( 6 ) 


can be rewritten as an integral over the plaquette variables U fiu (x) £ 517(3). The lattice Bianchi 
identity then automatically arises as the argument of a 5-function in the integrand of the partition 
function. It essentially constrains the product of plaquettes u covering a three-dimensional 
cube to the unit matrix in order to reduce the degrees of freedom (6 plaquettes per lattice site) 
to the original number (4 links per lattice site). On a four-dimensional lattice there are four 
types of three-dimensional cubes and thus four Bianchi identities. We want to emphasize that 
if a theory is expressed in terms of link variables the lattice Bianchi identities for the plaquette 
variables u are automatically fulfilled. In the case of an Abelian theory like 17(1) gauge theory 
the lattice Bianchi identity can be formulated also as a constraint for the plaquette phases. The 
sum of the phases 9^ £ (—7r,7r] of the six plaquette variables U /IU = e ie £ 17(1) measuring 
the flux out of a three-dimensional cube must be an integer multiple of 2ir. If the integer is 
non-zero, one speaks of a topological excitation called magnetic monopole [ 101. However, if one 
identifies the flux through a plaquette according to (|§) with sin 9 111 ], the total flux out of a 
three-dimensional cube is in general non-zero. The last statement is also true in the case of a 
non-Abelian theory without being in contradiction to the lattice Bianchi identity. We conclude 
that the sum of algebra elements on the right hand side of (|4]) is in general non-zero although 
the product of the corresponding group elements U is - according to the Bianchi identity - 
equal to the unit matrix. 

It is the purpose of this paper to study the influence of the magnetic current J m ,fi on the 
confining properties of the considered SU (3) lattice gauge theory. To be more precise, we define 
the following gauge invariant quantity 


I J, 


m,/x 


»l 



(7) 


and call it the length in colour space of the magnetic current defined in ([|). Large values 
of correspond to large fluctuations of the plaquettes covering the surface of a three- 

dimensional cube. We introduce a parameter A and modify the Wilson action ([l|) in the following 
way 

S = Sw + A^ | J m ,n{ x )\- (8) 

X,/-L 

The sum extends over all three-dimensional cubes on the lattice. The parameter A plays the role 
of a chemical potential for the magnetic current J m If A is chosen to be larger than zero, the 
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Figure 1: Expectation value of the magnetic current | J m \ ('9|) as a function of the inverse coupling 
/ 3 . The measurements were taken on a 12 4 -lattice. Errorbars are omitted since they are smaller 
than the symbols. 


magnetic current J m ^ will be suppressed. In ref. |T^] the action (|) with A = oo and with J m ^ 
being the 17(1) monopole currents |T(| was used to study four-dimensional 17(1) lattice gauge 
theory with monopoles removed. 

3 Results 

In order to study the influence of magnetic currents J m ^ on physical quantities we performed 
simulations using the action (j8[) with zero and non-zero chemical potential A. Calculations 
were done on symmetric lattices with size 8 4 and 12 4 and on asymmetric lattices with temporal 
extent N T = 4 and spatial extent of N a = 8 and 12. The case A = 0, i.e. the standard Wilson 
action, was simulated by using an overrelaxed pseudo heatbath algorithm. To generate gauge 
field configurations at non-zero chemical potential A we implemented a Metropolis algorithm. 
Because of the additional term in the action twelve three-dimensional cubes corresponding to 
colour magnetic currents J m ,n have to be taken into account besides the six plaquettes in the 
update of a single link. This procedure turned out to be very time consuming. We therefore 
decided to measure the physical quantities after every Metropolis sweep; this leads to large 
autocorrelation times especially in the critical regions of the finite temperature phase transition. 
As a consequence a large number of measurements up to 10 5 had to be taken into account in 
order to get a reliable estimate of the considered physical quantities. To estimate the errorbars 
of primary quantities such as Polyakov and Wilson loops the standard procedure of calculating 
the integrated autocorrelation time was used. For the error analysis of secondary quantities such 
as Creutz ratios a simple extension of the fundamental jackknife formula |l3j] was applied. 

We start with the presentation of results obtained on symmetric lattices of size 8 4 and 12 4 . 
Fig. a shows the expectation value of the magnetic current (j7[) averaged over the space-time 
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components and the four-dimensional lattice 


J rt 


1 

4 K 
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for a large range of /3-values and zero chemical potential A. One can see a typical crossover of 
(| J m \) between the strong and weak coupling region. For (3 going to infinity the magnetic current 
tends to zero displaying decreasing fluctuations of the six plaquettes covering a three-dimensional 
cube. In the strong coupling region the current (j|) surprisingly reaches its maximum at a finite 
value of (3 and slightly decreases for (3 going to zero. It is an interesting question whether the 
current density {\J m \} shows scaling behaviour which means (|J m |) oc a 3 (/3), where a{(3) is the 
lattice spacing given as a function of the bare coupling constant f3 which can be determined by 
the renormalization group equation. From the numerical data presented in fig. |l] it is clearly 
seen that there is no scaling behaviour of (\J m \). This may be explained by observing that 
there are only positive contributions to the current density (| J m \). Quantum fluctuations which 
are not of topological origin do not cancel. It should be emphazised that for the correlation 
function measuring the curl of the magnetic current in the vicinity of a static charge pair (see 
ref. U) this problem does not seem to occur. In this case, the magnetic current is projected 
to the Polyakov line which fixes a direction in colour space. Fluctuations corresponding to 
lattice artifacts contribute equally with positive and negative sign and thus cancel, whereas 
contributions of topological origin survive the averaging process. Unfortunately, the situation 
is different in the case of the operator \J m \ which measures the magnetic current density in 
the gluonic vacuum. We will come back to the question of scaling of this operator and its 
consequences to the continuum limit in the conclusions of this letter. 

To investigate the influence of the current J m)tl on physical quantities we calculated Wilson 
loops and Creutz ratios on a 8 4 -lattice at (3 = 6.0 and compared the results for zero and finite 
chemical potential A. In fig. || the ratio of Wilson loops obtained with A = 0.1 and A = 0.0 
(standard Wilson action) is shown. It is clearly seen that a suppression of magnetic currents 
manifests itself in a less disordered lattice: For all sizes of Wilson loops the considered ratio is 
larger than one and increases with increasing size of the loops. This means that the suppression 
of magnetic currents © corresponding to geometrical objects of the size of a unit cube especially 
influences the behaviour of objects being large in terms of the lattice spacing, i.e. large Wilson 
loops which determine the non-perturbative properties of the theory. To get a more quantitative 
insight we analyzed the behaviour of Creutz ratios x 


X(I, A) = - In 


W(I,I,X)W(I — 1,1 — 1, A) 
W(I — 1,1, X)W(I, I — 1, X) 


( 10 ) 


for the above discussed Wilson loops estimating the string tension. Since the ratio of Wilson 
loops shown in fig. |2] increases stronger than linearly, the ratio y(A = 0.1)/y(A = 0.0) decreases 
with the size of the loops. This behaviour is shown in fig. |j. The larger the size of Wilson loops 
the stronger is the influence of the suppression of magnetic currents on Creutz ratios. Taking 
the numerical results of \(4, A) as an asymptotic value for the string tension a in lattice units 
we are able to estimate the decreasing of the lattice spacing a 


a( X = 0.1) 


<t(A = 0.1) 
<t(A = 0.0) 


a(X = 0.0) 


0.79 o(A = 0.0). 


( 11 ) 
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Figure 2: Ratio of Wilson loops for chemical potentials A = 0.1 and A = 0.0. The ratio increases 
stronger than linearly with the size of the Wilson loops. 
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Figure 3: Ratio of Creutz ratios for chemical potentials A = 0.1 and A = 0.0 as a function of the 
linear extension of the considered Wilson loops. 
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We want to point out that for the chosen parameters f3 = 6.0 and A = 0.1 the average contribu¬ 
tion of the extra term A to the total action (||) is approximately 7%. 

In addition to zero temperature investigations we also performed numerical simulations at fi¬ 
nite temperature to study the influence of the magnetic current J mi#1 on the deconfinement phase 
transition. The order parameter of the finite temperature phase transition is the expectation 
value ( L) of the Polyakov loop 

N T 

L(x) = Tr ll U 4 (x,t). (12) 

t =l 

In the thermodynamic limit (L) is expected to be zero in the confined and non-zero in the 
deconhned phase. On a finite lattice, however, (L) equals zero for all inverse couplings (3 because 
of the non-zero probability of tunneling between different states related by the Z(3)-synnnetry. 
We thus measured the finite lattice “order parameter” (|L|) given by 

w> = dipE L (*)l>- < 13 ) 

JVfT X 

where the sum extends over the whole spatial lattice. We performed runs on lattices with size 
8 3 x 4 and 12 3 x 4 for chemical potentials A = 0.00 (standard Wilson action) and A = 0.05. The 
results for (|T|) as a function of the inverse coupling constant j3 are shown in fig. [| and ||. It is 
clearly seen that a non-zero chemical potential A shifts the phase transition towards a smaller 
critical value of (3. Hence, the addition of an extra term to the Wilson action in (|[) suppressing 
magnetic currents for A > 0 corresponds to a system with a larger effective inverse coupling 
constant. In other words, if configurations with large magnetic currents are removed from 
the path integral, the phase transition already occurs at smaller values of f3, and it is more 
favourable for the system to be in the Z(3 ) broken than in the Z( 3) symmetric phase. From 
fig. [| and |5] one can estimate the shift of the phase transition. For A = 0.05 we find that A (3 
is approximately —0.2, whereas the current contribution A^ 1(1 | to the total action (|j) is 
about 4% for (3 £ [5.3, 5.7]. 

Besides the Polyakov loop we also measured the distribution of the real part of the plaquette 
trace. In fig. |] we show a comparison of plaquette distributions on the 12 3 x 4-lattice at (3 = 5.6 
for A = 0.00 and A = 0.05. The corresponding numerical values of (\L\) are displayed in fig. J5| 
As expected the suppression of magnetic currents causes the plaquettes to be closer to the unit 
matrix. Plaquette distributions corresponding to the same value of (|L|) but different values of 
[3 for A = 0.00 and A = 0.05 cannot be distinguished. 

By choosing the chemical potential A large enough the confinement phase can be completely 
removed for positive values of (3. Starting from an ordered lattice and choosing A = 10.0 the 
system is not able to evolve from the deconfined into the disordered confined phase, also not for 
(3 = 0.0, see fig. [?]. This shows that for the considered lattice sizes 8 3 x 4 and 12 3 x 4 there is no 
confined phase without presence of magnetic currents for positive values of (3. 

4 Conclusion and discussion 

In this letter we studied the influence of colour magnetic currents @ originally defined in ref. || 
on the confining properties of SU(3) lattice gauge theory. We introduced a chemical potential 
A and added an extra term to the standard Wilson action suppressing large magnetic currents. 
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Figure 4: The expectation value of the Polyakov loop as a function of the inverse coupling (3 
for chemical potentials A = 0.00 and A = 0.05 for a 8 3 x 4-lattice. In the case of A ^ 0 the 
critical /3 is shifted towards a smaller value. The symbols (diamonds for A = 0.00 and triangles 
for A = 0.05) denote the numerical results of Monte Carlo runs at different values of (3 whereas 
the continuous curves were computed with a multi-histogram analysis [14, 151. 



Figure 5: The same as in fig. [| for a 12 3 x 4-lattice except that the curves are linear interpolations 
between the results of Monte Carlo runs drawn to guide the eye. 
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Figure 6: Distributions of the real part of the plaquette trace for (3 = 5.6 and the A- values 0.00 
(dashed line) and 0.05 (solid line). According to fig. [5| the lattice is in the confined phase for 
A = 0.00 whereas it is in the deconfined phase for A = 0.05. 
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Figure 7: The expectation value of the Polyakov loop as a function of the inverse coupling j3 for 
chemical potential A = 10.0 for a 8 3 x 4-lattice (left plot) and a 12 3 x 4-lattice (right plot). At 
this value of A there is no confined phase left for positive values of the inverse coupling (3. 
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With this modified action we performed numerical simulations at zero and finite temperature 
lattices. We found that a partial absence of magnetic currents leads to a drastic decrease in 
the confining string tension. At finite temperature the suppression of magnetic currents shifts 
the phase transition towards smaller values of the inverse coupling constant /3. Additionally, 
for the considered lattice sizes it was shown that there exists a finite chemical potential A for 
which the confined phase is completely removed for positive values of (5. We conclude that the 
magnetic currents strongly influence the confining properties of our considered SU{ 3) lattice 
gauge theory. However, an interpretation of these currents in terms of physical objects is not 
directly possible since we do not observe scaling behaviour of the current density (|9j) . As already 
mentioned above, the reason for this might be that there are only positive contributions to the 
expectation value (| J m \) and that fluctuations which are not of topological origin do not cancel. 
According to our opinion, on a lattice with finite spacing there are in general two different types 
of contributions to the magnetic current One type of contributions is due to the finite 

lattice spacing. They can be regarded as lattice artifacts and vanish in the continuum. Hence, 
they will not contribute to the continuum limit of the expectation value (| J m \). The other type 
of contributions is of topological origin. They belong to gauge field configurations which become 
singular in particular gauges in the continuum and thus could guarantee that the operator J m 
does not become trivial in the continuum limit. One candidate for such configurations is for 
instance the Euclidean pure gauge theory analogue of the ’t Hooft-Polyakov-like monopole in 
the BPS limit. That this monopole configuration could play an important role in the dual 
superconductor picture of confinement was discussed in ref. |l^]. Moreover, in ref. [17| it was 
argued that in an appropriate gauge such a monopole configuration becomes singular. Analogous 
arguments should be valid for another type of gauge field configurations, magnetic vortices. 
These were suggested to play a crucial role in the confinement mechanism [18|, an idea which 
gained support from recent lattice calculations m- 

To summarize, we believe that the topological contributions to the magnetic currents strongly 
influence the confining properties of QCD. Unfortunately, the magnetic current operator is not 
able to identify its topological sources which determine the behaviour of its expectation value 
in the continuum. Thus, we can only speculate about this question. But for cubes of finite size, 
we know that the topological objects lead to magnetic currents behaving as predicted by the 
dual superconductor picture. 
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